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I. INTRODUCTION 


Spontaneous spin-up or flow generation is of great interest in many fields of scientific 
research, such as the jet-stream formation in the atmosphere [l|, the giant zonal belt on the 
Jupiter the flow pattern generated in a turbulent neutral fluid jfj, the dynamo effect 




m 


a conducting fluid |4|, the particle acceleration in solar flares jf|, and the zonal flows in 
a tokamak magnetic fusion plasma jg]. The toroidal rotation is crucially important for the 
tokamak fusion reactor, a magnetic fusion torus, such as the International Thermonuclear 

n nn 

Experimental Reactor (ITER) [7], since it can stabilize the macroscopic instabilities [8|, 9], 
and regulate the microturbulence fiol 13 ]. Recent tokamak experiments suggest that an 


intrinsic toroidal rotation exists 


20], due to an anomalous momentum source in the edge 


region, which can not be explained by the effects of momentum pinch 


21 


22]; the effects of 


the residual Reynolds stress due to parallel symmetry-breaking on the parallel momentum 
generation 23] and the effects of toroidal symmetry-breaking on the toroidal acceleration 
24 ] have been discussed. 

In this paper, we shall show that the intrinsic toroidal rotation can be generated by the 
turbulent (anomalous) radial particle flux, which is reminiscent of the bootstrap current 
generated by the neoclassical particle radial diffusion 25|, |26j; the proposed mechanism of 
toroidal momentum generation can be explained by the electron-ion asymmetric anomalous 
viscous damping, which does not need the spatial symmetry-breaking, such as the toroidal 
symmetry-breaking or the parallel symmetry-breaking. This paper is organized as follows. 
In Sec. II, we review the fundamental theory and the related experiments; in Sec. Ill, we 
present the theoretical analysis; in Sec. IV, we discuss the applications of the proposed 
model; in Sec. V, we present the summary and discussions. 


II. REVIEW OF THE FUNDAMENTALS 

II. 1. Fundamental theory 

The equilibrium magnetic field of a tokamak is written as B = I (r) V£ + V£ x (r), 
with / = RBt, Bt the toroidal magnetic field, R the major radius. The poloidal magnetic 
flux is ^ (r), with r essentially the minor radius of the torus. In the large-aspect-ratio limit 
(r <C R), the poloidal magnetic field is given by Bp = i^'/R, with the prime denoting the 
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derivative with respect to r. ( is the toroidal angle. 

The two-fluid mean-field momentum equation for a simple electron-deuteron plasma can 
be written as [27], 


d t ( n s m s u a ) + V-(n s m s u s u s + p s I + U s n + n SiQ ) = n s e s [-V<f> (0) + E ( V( + u s x B] + F s , 

( 1 ) 

where the subscript s = e,i denotes the charged particle species, m and e are particle mass 
and charge, respectively, n and u are density and fluid velocity, respectively. Note that the 
quasi-neutrality condition demands that n e = rii = n. —V$(0) and E0V ( = E T are the 
radial electric field and the toroidal inductive electric field, respectively, p is the pressure. 
n„ is the CGL form of the stress tensor, which includes the neoclassical effects [27J. II a is 
the anomalous viscosity tensor, which accounts for the anomalous transport of momentum 


due to electrostatic turbulence, therefore, it includes the usual Reynolds stress term 24], 


28]. 


F is the collisional friction force. In Eq. (JT[) , we have neglected the source term due to the 
external momentum injection. Note that all quantities in Eq. (0Q) are turbulence-ensemble 
averaged. 

Acting i? 2 V<0 on Eq. (JT]) and averaging the resulting equation over the magnetic flux 
surface of a tokamak yields the toroidal angular momentum equation 27] 


d t {n s m s u StC ) + (V ( n s m s u r s u sX + fl^)) - ((F S)C ) + n s e s E ( ) = e s Y r s (2) 


dr 


where V (r) is the volume enclosed in the magnetic flux surface labeled by r. Note that 
he CGL form of the neoclassical viscosity tensor does not survive the flux surface average 


27]. The radial particle flux is T r = n (r) ( u r ). Note that u r = Vr • u, = A 2 VC ■ u, and 


n£ = Vr • 7T • i? 2 VC- 

A simple interpretation of Eq. ([2]) is as follows. The radial particle flux can be decomposed 
into two components: the neoclassical radial particle flux T n and the anomalous radial 
particle flux T a . 


r r = r r + r r 

.«? fit n. 1 fit 


(3) 


The neoclassical radial particle flux is driven by the toroidal collisional friction force and 
the toroidal inductive electric field, 


e s r r J)' = - ((F a>< ) + n s e s E ( ); 


(4) 
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the anomalous radial particle flux is driven by the anomalous toroidal force written in terms 
of the anomalous toroidal viscosity and the inertial term, 


1 d 


e s r3,00' = d t ( n s m s u s , ( ) + — — (V' ( n s m s u r s u sX + II^)) 


(5) 


This interpretation of the anomalous radial par 


theories (Sec. 4 of Ref. 


28 


and Sec. II. C of Ref. 


icle flux is consistent with the previous 


24|); the toroidal momentum equation used 


in Ref. [24j includes an additional term accounting for the toroidal turbulence acceleration, 
which depends on the toroidal symmetry breaking and shall be ignored in this paper. Fol- 

, one concludes that although Eq. (J2]) was originally derived for an 
is also applicable for a turbulent plasma with the effects of Reynolds 


lowing Ref. 


24 and 28 


H it i 


quiescent plasma 
stress included in the anomalous viscosity term. 

It is useful to make some comments on the neoclassical bootstrap current in the large 


aspect-ratio limit. The collisional friction force is given by 


29] 


Fi F ie F e i U e iTiTTl e (Ui ’U'e) ? 


( 6 ) 


with v ei the electron-ion collision frequency. 

The neoclassical radial particle radial flux, which satisfies the auto-ambipolarity condi¬ 
tion, is given by [27] 

dn 
dr 


r r = T r ~ 

i,n e,n 


-ftVeiP, 


ep- 


ftnE^'/R 2 , 


(7) 


with f t the trapping ratio. p eP = \/2m e T/e e B P is the electron poloidal gyro-radius, with T 
the temperature of the plasma. 

Substituting Eq. (|6]), Eq. (j7J) into Eq. (J4]), one finds 


(n e ei (Ui 



~ f ‘l P d i + 


( 8 ) 


with a s the Spitzer conductivity. This shows the well-known bootstrap current and the 
neoclassical reduction of the Spitzer conductivity. 


II.2. Experiments related to the friction term 

Both of the toroidal rotation and the toroidal current are crucially important for a toka- 
mak fusion reactor. Fortunately, the toroidal current can be partly provided by the bootstrap 
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current driven by the radial neoclassical particle flux. The neoclassical bootstrap current 
and electric conductivity depends on the friction term, which has been discussed in the 
above. The friction term is important in the toroidal momentum equation. Therefore, it is 
useful to examine the experiments related to the friction term. 

Although generally tokamak experiments indicate that the radial particle flux is anoma¬ 
lous (typically 30 times larger than the friction-generated neoclassical radial particle flux), 
one can attribute the anomaly of the radial particle flux to an anomalous toroidal force (due 
to anomalous toroidal viscosity). However, it should be pointed out that many experiments 
have been carried out to test the neoclassical theory of bootstrap current and the Spitzer 
conductivity [30H33]. These experimental observations indicate that the toroidal current is 
related to the toroidal inductive electric field and the radial pressure gradient in the way 
predicted by the neoclassical theory based on the collisional friction. Although there are 
theoretical arguments that the bootstrap current may be enhanced by the effect of radial 


electric field 


34!] . or by the effect of turbulence scattering 35j, the order of magnitude of the 


bootstrap current and hence the friction does not change. Therefore, one concludes from 
the experimental observations that the toroidal friction term is not anomalous; the toroidal 
friction force acted on one species is balanced by the poloidal magnetic field crossed with the 
radial current carried by its neoclassical flux. Eq. (J8]) has been validated by experiments. 


II.3. Experiments related to the intrinsic rotation 


Intrinsic toroidal rotation has been observed in many tokamak experiments 14H201], which 


indicate that significant co-current toroidal rotation exists in tokamaks without any external 
momentum injection; this important phenomena is difficult to understand, especially when 
considering that a tokamak is essentially a toroidal symmetric system. 

In C-mod, during a transient phase of density ramp-up, a toroidal counter-current accel¬ 
eration is observed for both the ion-cyclotron-resonance heated plasma and the lower-hybrid 
wave-driven plasma, and a co-current acceleration is observed during a density ramp-down 
phase 18]. The most convincing evidence of the anomalous torque and hence the sponta¬ 


neous rotation may be the recent measurement of the anomalous co-current torque reported 
by the DIII-D team ill 2 o|; the anomalous torque density is in the co-current direction, 
with its peak value ~ 0.7Nm/m 3 located at the edge. 
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The observation of an anomalous toroidal torque in the edge o: 
with the observation of spontaneous toroidal rotation in C-mode 
is used to interpret the toroidal rotation. 


DIII-D 


3,20] 


is consistent 


16], where a model equation 


dtP H— d r 


r [ ~X(drP - -V C P 


= 0 , 


P(r = a) = P 0 , 


(9a) 

(9b) 


with P the toroidal momentum of the plasma, a the minor radius of the device, and v c the 
momentum pinch velocity. Note that a finite toroidal momentum at the edge in Eq. (19b|) 
implies a source at the edge when a no-slip boundary condition is applied. 


III. THEORETICAL ANALYSIS 


Consider a large-aspect-ratio tokamak. Eq. ([2]) is reduced to 


1 d 


d t ( n s m s U s ) H-- < rn s m s 

r dr 


d 


K + «U) u, - x„( T u, 


- R (F S £ + n s e s E() = e s Y r s B P , 

( 10 ) 

where U s = RV( ■ u s is the toroidal velocity. Xs,c the momentum diffusivity. u r s mp is the 
radial pinch velocity of the toroidal momentum. Note that the anomalous toroidal viscosity 


is written in a momentum diffusion term anc 
The ion momentum pinch is given by 


a momentum pinch term. 


«„m P = 


( 11 ) 


with a = 2, 4, to take into account of the momentum pinch due to the Coriolis force 21 
For simplicity, we shall assume that this toroidal momentum pinch-diffusion relation is also 
applied to electrons. 

The neoclassical particle flux and the anomalous particle flux are written as 

1 


pr 

s ’ n " e s D 


s^jp 


-R ( P.s.C, + Tl's&sPc,') , 


-pr _ 

~ e P B 


1 „ . TT x lid 

d t ( n s m s U s ) + -pT-— 1 rn s m s 


{w,-%x„c)u,- x , c d-u. 


i s Dp e s Bprdr 

The ion radial particle transport is schematically shown in Fig. 1. 


( 12 ) 

(13) 
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The ambipolarity condition, which is derived from Ampere’s law, demands that 


e e Tf + eW r = 0. 


(14) 


Since the neoclassical particle flux is auto-ambipolar, the ambipolarity condition is reduced 
to 


r 


r 

e,a 



(15) 


For typical tokamak experimental parameters, U,. ~ 3 x 10 4 m/s, U e ~ 10 6 m/s, 


rrijUi 

m e U e 


10 2 . 


Following Eq. (US]) and Eq. (1T5|) . one concludes that 


(16) 


^ ~ 10 2 , 

Xi, c 


which indicates an electron-ion asymmetry of anomalous toroidal viscous damping. 


(17) 


III.l. Relaxation process 


For typical tokamak experimental parameters, v ei ~ 10 5 /s; the ion toroidal momentum 
confinement time is ~ a 2 / Xi, C rs "'’ , 10 x s; T r /na = u r /a ~ 1 m/s. Following Eq. (117j) . one 
finds T e £ a 2 /Xe, 

C rs ~' 10 3 s. Using these parameters, one finds that the electron momentum 
convection term due to the particle flux can be safely dropped. 

We write down the electron toroidal momentum equation 

£4 j-rie [Vtim e (U, : - Ui) + e e E T ] = e e (r'„ + C J B P . 

(18) 

We observe that 

1 /Vei < Te,C < T i,Ci ( 19 ) 

which indicates that the collisional relaxation of U e is faster than the turbulence relaxation. 

On the short time scale of l/r e j, U e is quickly relaxed through the neoclassical collisional 
process, 

v ei n e m e (U e - Ui) + n e e e E T = e e T r e n B P , (20) 

which is consistent with Eq. (jHJ) and Eq. (1T2]) . 


1 d 


d t ( n e m e U e )+rn e m e Xe£ 


a 

R 


d_ 

dr 






On the longer time scale of r eX , U e is relaxed through the turbulent process to generate 
the anomalous radial electron flux 


1 d 

r dr 


rn e m e 


a tt d TT 

R Xe, C Ue XeX^e 


= eeK.aBp, 


( 21 ) 


On the time scale of T t ^, the ion toroidal momentum is relaxed to generate the anomalous 
radial ion flux 


d t (riimiUi) + rriirrii 

r dr 


a TT d TT 

R Xi,cUi x« d U< 


= or i a B P . 


( 22 ) 


III.2. Single-fluid toroidal momentum equation 


The single-fluid toroidal velocity V is defined by pV = nMV = n(mi + m e )V = n(rriiUi + 
m e U e ). The single-fluid toroidal momentum equation can be found by adding the electron- 
version of Eq. ([TOP to its ion-version and using the ambipolarity condition. Since Eq. (HO 
can be written as Eq. (TT3l) and Eq. ([T2lh and the ambipolarity condition is reduced to Eq. 

one can write down the single-fluid equation through combining Eq. (fl^jl . Eq. (TTTT1) 
and Eq. (TT2|) . 

1 d 


dt(pV)+-^- \ rp 
r dr 


a \ d 

U ' ~ R Xi ’ C ' ~ X *’ C 


dr 


V =- 


r dr 


mm e { Xe ,c ~ Xi,c ) ( 

' (23)' 


which is the exact result of the ambipolarity condition. 

A single-fluid toroidal momentum equation consistent with Eq. (I9ap . which has been 
widely used in previous works |16|, 119|, [20] on the transport of toroidal momentum, can be 
found from Eq. (1231) by simply dropping its right-hand side. If X e,( = Xe£, this would have 
been true. If Xe,£ ~ Xe,o Eq. (JHIJ would have been a fairly good approximation for typical 
experiments for which m,iUi m e U e . However, according to the above discussions, for 
typical experiments, there this a strong electron-ion asymmetry of toroidal viscus damping 
[see, Eq. (JTTD]. 

Using Eq. (Eli and Eq. (j2B . one finds that Eq. (j23l) can be written as 

d 


dt ( pV) + < r 

r dr 


(v - pV - p Xu -V 

Using Eq. (fT5p . and defining the anomaly factor 


—e T r B 


(24) 


pr _ f ur 

1 a — J a L ? 


(25) 
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one finds the single-fluid toroidal momentum equation 

1 d 


dt ( pV ) + < r 

r dr 


(“' - v«) pV ~ mx i v 


= U,B p r 


(26) 


Note that typical experimental parameters indicate that f a < 1. 

General remarks should be made on Eq. (J26]). The right-hand side identifies the anoma¬ 
lous momentum source observed in experiments. Note that usually the particle source is 
in the edge of a tokamak due to the shallow penetration of the fueling; therefore, in a 
steady-state, the particle flux is non-zero only in the edge region. This explains why the 


experimental observations indica 
is located in the edge region 


;e that usually an anomalous toroidal momentum source 
17|. The source term, the right-hand side of Eq. (|24|) 
represents the intrinsic toroidal rotation driven by the anomalous radial particle flux, which 
is different from the previous models which need a spatial (essentially toroidal) symmetry¬ 


breaking effect 


23 


24] 


The physical mechanism of spontaneous momentum generation behind Eq. (126|) can be 
understood by examining Eq. (fTO . The friction forces acted on electrons and ions exactly 
cancel out due to momentum conservation, and the forces due to the toroidal inductive 
electric field cancel out through the quasi-neutrality condition. Although the toroidal com¬ 
ponent of magnetic Lorentz forces ( e s T r s Bp ) cancel out due to the ambipolarity condition, 
the electron toroidal viscous damping rate (Xe,() I s larger than the ion toroidal viscous damp¬ 
ing rate (Xi,c) [ see ; Eq. (117]) ]. this electron-ion asymmetry of viscous damping can generate 
a net momentum for the two-component system. 

The spontaneous ion toroidal rotation can also be understood in a two-fluid picture. 
From the ion-version of Eq. (fT3]h one can see that the anomalous toroidal ion viscous force 
is balanced by the product of the anomalous radial ion current and the poloidal magnetic 
field, which is schematically shown in Fig. 1. 


IV. APPLICATIONS 


If one assumes a constant value of Xi, C, the analytic solution of Eq. (1261) with a no-slip 
boundary condition can be readily found, and the toroidal velocity is given by 


U (r) = [c (a) — c (r)] exp 


v r + v r 
7 u* ' U'z.mp 

dr --— 

x« 


(27) 
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where a is the minor radius of the boundary, 


c(r) = / dr- 


exp — / dr 


Iq rnrriiXix 


S ( (r) = / rdrneiif B 


n r + ul mp 

Xi,c 


p- 

Jo 

u r can be determined as follows. 

Particle balance demands 

2nmu r = / 2irrdrSog (r), 

Jo 

with Sog (r) the particle source and So a constant. So can be determined by 

[nu r ] a 2na = 7ia 2 n/T p , 


(28) 

(29) 


(30) 


(31) 


with n the volume averaged density. 

Typical parameters of an ITER-like tokamak are as follows. R/a = 6.2m/2m, Bt = 5.3T, 
e* = 1.6 x io- 19 g, rrii = 2.5 m p . The energy confinement time is te ~ 3.5s, the particle 
confinement time is r p ~ 10s, which is roughly 3 times of the energy confinement time. We 


shall assume n = n 0 


0.5 


1 — (^) + 0.1n 0 , with n 0 = 1.5 x 10 2 °/m 3 , and q (r) = = 


1.1 


1 + 2.5 (p' 


Since the toroidal momentum confinement time is usually similar to the 


energy confinement time 


0 . 


we shall assume Xi£ 


5 . 8 te ‘ 


A Gaussian type of particle 


source profile is assumed, which is non-zero in the edge region 0.9 < r/a < 1. 

The particle source Sog (r) and the toroidal force density neiU r Bp are shown in Fig. 2. 

With the above parameters, the toroidal rotation speed for different momentum diffusiv- 
ities (x^) with different momentum pinch factors ( a ) are shown in Fig. 3. 

Note that the solution of the core rotation (r < 0.9) shown in Fig. 3 can also be 
obtained by solving the source-free version of Eq. (1261) and adjusting the prescribed boundary 
condition V (r — 0.9) = Uq, which is similar to the assumption used in Ref. [16 


Fig. 3 clearly demonstrates that a significant spontaneous toroidal rotation can be driven 
by the anomalous particle diffusion in an ITER-like tokamak fusion reactor, which is typically 
from 200 km/s in the edge to 500 km/s in the core, and this bootstrap rotation may be 
sufficient to stabilize the resistive wall mode [l7, 36]: it is interesting to note that the 
spontaneous co-current toroidal rotation predicted by this model roughly agrees with the 
value extrapolated from the present experimental data 


0 


17]. The effect of momentum pinch 


, 22] can be seen by contrasting Fig. 3 (a, b) to Fig. 3(c). 
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It is of interest to compare the proposed model with the recent experimental observations 
of the spontaneous rotation. In a steady-state plasma fueled at the edge, the particle flux is 
radially outward, thus produces a co-current torque at the edge; this prediction qualitatively 
agrees with the experimental observations of the spontaneous rotation of the H-mode plasma 
LL9[. The direction of the spontaneous rotation of the L-mode plasma is complicated 
17], and the counter-current spontaneous rotation of the L-mode plasma may be due to the 


neoclassical effects 


20 


37]. During a transient phase of density ramping-up, a radially inward 


particle flux is expected, which should generate a counter-current torque, and this prediction 
qualitatively agrees with the experimental observations of the coimter-current acceleration 
during the density ramping-up phase of both the ion-cyclotron-resonance heated plasma and 
the lower-hybrid wave-driven plasma 18]. The most convincing evidence of the anomalous 
torque and hence the spontaneous rotation may be the recent measurement of the anomalous 
co-current torque reported by the DIII-D team juj. 20]; the anomalous torque density is in 
the co-current direction, with its peak value ~ 0.7Nm/m 3 located at the edge. The main 
parameters of the DIII-D tokamak is R/a = 1.6m/0.6m. For the experiments discussed 


in Ref. 


19 


one may assume that the volume-averaged electron density is ~ 3 x 10 19 /m 3 , 


the particle confinement time is r p ~ 0.3s, and the poloidal magnetic field at the edge is 
~ 0.15T. With these parameters, following the above analysis, one finds that the peak value 
of the torque density at the edge is ~ 1 Nm/m 3 , which agrees well with the experimental 


value 


19 


20 ]. 


V. SUMMARY AND DISCUSSIONS 


In summary, we have proposed a theoretical model of spontaneous toroidal rotation driven 
by the turbulent radial particle flux through the electron-ion asymmetric viscous damping. 
The proposed model explains the anomalous co-current momentum source observed in the 
tokamak edge, and predicts a significant toroidal rotation in ITER, which may mitigate 
the needs of momentum injection and modify the design of a tokamak fusion reactor. It 
is interesting to note that both the toroidal current and the toroidal rotation are crucially 
important for the steady-state operation of tokamak fusion reactor 361]; the bootstrap cur¬ 
rent is driven by the neoclassical particle diffusion and the intrinsic toroidal ion rotation is 
driven by the anomalous radial particle flux. 
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The mechanism of spontaneous momentum generation in a two-component system with 
the asymmetric damping can be simply understood as follows. Consider a two-component 
system given by 


d t P\ — —V\P i — i / 0 (-Pi — P2) + P-i (32a) 

d t P 2 = -u 2 P 2 - v 0 (P 2 - Pi) - F, (32b) 

where P s is the momentum of each component; v s is the viscous damping rate; vq is the inter¬ 
component friction rate; F is the driving force acted on component 1 , which is opposite to 
the driving of component 2. The initial values are set as P\ = 0 = P2. The final steady-state 
total momentum of the system is given by 

Py + ft = — FPr 1 F - (33) 

viv 2 + (z/i + i/ 2 ) i/o 

which clearly demonstrates a net momentum generated by the asymmetry of viscous damp¬ 
ing, when the driving forces cancel out. 

The physical mechanism of the electron-ion asymmetry of toroidal viscous damping is 
still an open issue, although it appears to be consistent with experimental observations. 
This paper proposes a new way to look into the physics of spontaneous toroidal spin-up; 
the present researches concentrate on the ion viscosity and the effects of possible parallel 
or toroidal symmetry breaking on the ion dynamics, this paper suggests that the electron 
viscous damping is important. A possible explanation of the electron-ion asymmetric toroidal 
viscous damping is as follows. Consider the trapped electron mode. Since the mode is 
resonant with electrons instead of ions, there should be a resonant electron momentum 
diffusion and no resonant diffusion of the ion momentum; this shall be left for the future 
researches. 

Figure Captions: 

Fig. 1 . Ion radial flux driven by anomalous toroidal viscous force and the neoclassical 
toroidal friction. 

Fig. 2. Particle source and toroidal force density. 

Fig. 3. Toroidal rotation speed for different momentum diffusivities. Xo — ci 2 /5.8te ■ (a) 
a = 4; (b) a = 2; (c) a = 0. 
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